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Abstract 


Results of a search for new phenomena in events with at least three photons are reported. 
Data from proton-proton collisions at a centre-of-mass energy of 8 TeV, corresponding to an 
integrated luminosity of 20.3 fb”', were collected with the ATLAS detector at the LHC. The 
observed data are well described by the Standard Model. Limits at the 95% confidence level 
on new phenomena are presented based on the rate of events in an inclusive signal region 
and a restricted signal region targeting the rare decay Z ^ 3y, as well as di-photon and 
tri-photon resonance searches. For a Standard Model-like Higgs boson with a mass of 125 
GeV decaying to four photons via a pair of intermediate pseudoscalar particles (a), limits are 
found to be cr X BR(/i ^ aa) X BR(a ^ yy)^ < IO'^ctsm for 10 GeV < < 62 GeV. Limits 

are also presented for Higgs boson-like scalars (H) for mh > 125 GeV, and for a Z' decaying 
to three photons via Z' a + y ^ 3y. Additionally, the observed limit on the branching 
ratio of the Z boson decay to three photons is found to be BR(Z —> 3y) < 2.2 X 10”®, a result 
five times stronger than the previous result from LER 
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1 Introduction 


Many extensions of the Standard Model (SM) include phenomena that can result in final states consisting 
of three or more photons. Extensions of the SM scalar sector [1-5], for example, often include pseudo¬ 
scalar particles (a) with couplings to the Higgs boson [6,7] (h) and branching ratios into photons that 
would be visible at the LHC, in addition to scalars (H) with masses different from the SM-like Higgs 
boson of mh = 125 GeV that can also decay via H ^ aa ^ 4y. Other models feature additional vector 
gauge bosons that can decay to a photon and a new pseudoscalar boson, a, with the subsequent decay of 
the a into a pair of photons, resulting in a three-photon final state [8]. Moreover, in the SM, the Z boson 
can decay to three photons via a loop of W- bosons or fermions. The decay is heavily suppressed and the 
branching ratio is predicted to be ~5xl0~'° [9]. The current most stringent bound on this process comes 
from the L3 Collaboration, which placed a limit of BR(Z ^ 3y) < 10“^ [10]. The ATLAS detector has 
collected ~10^ Z boson events, and thus an observation of this decay would indicate an enhancement of 
this decay rate and could be evidence of phenomena not predicted by the SM. Leynman diagrams for 
some of these beyond-the-Standard Model (BSM) and rare SM scenarios are shown in Ligure 1. 



Figure 1: Feynman diagrams for possible beyond-the-Standard Model (top) and rare Standard Model (bottom) scen¬ 
arios that result in final states with at least three photons. 


To ensure sensitivity to these and other possible rare SM and BSM scenarios, an inclusive three-photon 
search is performed using 20.3 fb“^ of LHC proton-proton collisions collected by the ATLAS detector 
in 2012 at a centre-of-mass energy of 8 TeV. Such a model-independent search is the first of its kind, 
as are the interpretations for a Higgs boson decaying to four photons via two intermediate pseudoscalar 
a particles (for a Higgs boson of m/, = 125 GeV and for Higgs-like scalars of higher masses) and for 
three-photon resonances corresponding to a new vector gauge boson. 

The dominant backgrounds include the irreducible component with three or more prompt photons, as 
well as the reducible components consisting of combinations of photons and electrons or hadronic jets 
misidentified as photons. The contributions from events with jets which are misidentified as photons 
are calculated from data-driven methods, while simulation is used to estimate the contributions from the 


2 



irreducible background and the reducible background originating from electroweak processes that lead to 
electrons which are misidentified as photons in the detector. Collision data is used to derive corrections 
to the probability obtained from simulation that electrons are misidentified as photons. 


2 The ATLAS detector 

The ATLAS experiment [11] at the LHC is a multi-purpose particle detector with a forward-backward 
symmetric cylindrical geometry and a near An coverage in solid angle. ' It consists of an inner tracking 
detector surrounded by a thin superconducting solenoid providing a 2 T axial magnetic field, elecfromag- 
nefic (EM) and hadronic calorimefers, and a muon specfromefer. The inner fracking defecfor covers fhe 
pseudorapidify range \r]\ < 2.5. If consisfs of silicon pixel, silicon micro-strip, and transition radiation 
tracking detectors. Lead/liquid-argon (LAr) sampling calorimeters provide EM energy measurements 
with high granularity. A hadronic (iron/scintillator-tile) calorimeter covers the central pseudorapidity 
range (I77I < 1.7). The end-cap and forward regions are instrumented with EAr calorimeters for both EM 
and hadronic energy measurements up to I77I = 4.9. The muon spectrometer surrounds the calorimeters 
and is based on three large air-core toroid superconducting magnets with eight coils each. Its bending 
power ranges from 2.0 to 7.5 T m. It includes a system of precision tracking chambers and fast detectors 
for triggering. A three-level trigger system is used to select events. The first-level trigger is implemented 
in hardware and uses a subset of the detector information to reduce the accepted rate to at most 75 kHz. 
This is followed by two software-based trigger levels that together reduce the accepted event rate to 
400 Hz on average depending on the data-taking conditions during 2012. 


3 Event and object selection 

This search utilises a three-photon trigger that places a minimum requirement on the photon momentum 
in the plane transverse to the beam axis (transverse momentum, or pt) of 15 GeV, applied on three photon 
candidates in the EM calorimeter. Each candidate is additionally required to satisfy a set of loose photon 
identification criteria [12]. Stringent detector and data quality criteria are applied offline. Evenfs are 
required fo confain af leasf one inferacfion verfex, wifh no additional vertex requiremenfs. 

Photon candidafes musf safisfy a pseudorapidify requiremenf of Ipl < 2.37, excluding fhe fransifion region 
befween fhe barrel and end-cap of 1.37 < \r]\ < 1.52, and musf safisfy requiremenfs on fhe shape of fhe 
energy deposif in fhe calorimefer. A phofon candidafe is rejecfed if fhe barycenfre of ifs energy deposif 
is wifhin a cone of A/? = '\/(A77)2"+7A0)2 < 0.15 around fhe barycenfre of fhe energy deposif of a higher 
Pt photon candidafe. Einally, selecfed phofon candidafes are required fo safisfy a more sfringenf sef 
of identification criferia, known as tight [12]. Phofon isolafion is defined by fhe amounf of fransverse 
energy, deposifed in fhe EM calorimefer wifhin a cone of size A/? around fhe phofon candidafe, 
excluding fhe energy of fhe phofon candidafe ifself. If is a powerful means of disfinguishing befween 
phofons and hadronic jefs misidenfified as phofons, since fhe energy clusfers deposifed by phofons in fhe 
EM calorimefer fend fo be narrower in fhe fransverse direcfion fhan fhose deposifed by jefs. Because 


* ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in the centre of the detector 
and the z-axis along the beam pipe. The x-axis points from the IP to the centre of the LHC ring, and the ^-axis points 
upwards. Cylindrical coordinates (r, 0) are used in the transverse plane, (f being the azimuthal angle around the beam pipe. 
The pseudorapidity is defined in terms of the polar angle 0 as ?/ = - In tan(0/2). 
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minimum-bias proton-proton interactions in the same or nearby bunch crossings (pileup) can affect the 
calculated photon isolation energy, a correction is applied based on an event-by-event energy density 
pileup estimation. This search uses an isolation cone of size AR < 0.4, and a correction to the value 
of a photon candidate is made when another photon candidate passing the tight identification criteria is 
found within an annulus of 0.15 < AR < 0.4 around the photon candidate. The correction consists of 
subtracting the pj value of the other photon candidate found within the annulus from the value of the 
photon candidate under consideration. The final isolation criterion is Ej° < 4 GeV. 

Events in the inclusive signal region are required to have at least three tightly identified and isolated 
photon candidates, where the two photon candidates with the highest transverse momentum must have 
Pj > 22 GeV while the third highest must have pj > 17 GeV. The restricted signal region, targeting 
the rare decay Z —> 3y, is a subset of these events where an additional criterion of 80 GeV < m^y < 
100 GeV is placed on the invariant mass of the three-photon system. The signal regions are supplemented 
by several control regions, where at least one of the photon candidates fails the isolation requirement. 


4 Simulated event samples 

Simulated event samples are used to estimate several SM background processes in the search for excesses 
in the inclusive signal regions, as well as to model signal predictions for both the inclusive searches 
and the resonance searches for the specific BSM scenarios considered here. The simulated BSM sig¬ 
nal samples are also used to define a fiducial region for which the search criteria are largely model- 
independent. 

4.1 Simulated backgrounds 

The SM two-photon process is an irreducible background to a three-photon search because a third photon 
may arise from minimum-bias proton-proton interactions in the same bunch crossing. The SM two- 
photon background is simulated with Pythia 8 [13], and the three- and four-photon backgrounds are 
simulated with MadGraph 5 [14], with Pythia 8 used for fragmentation and hadronisation. The produc¬ 
tion of two, three and four photons in the SM contains large contributions from higher-order Feynman 
diagrams. Thus, the two-, three- and four-photon simulated event samples calculated at leading order 
(LO) are multiplied by factors (“A'-factors") determined from studies with generators that include next- 
to-leading order (NLO) contributions, namely MCFM 6.8 [15] and VBFNFO 2.7.0 [16-18], using the 
parton distribution function (PDF) sets CTEQ6F1 [19] for the FO cross sections and CTIO [20] and 
MSTW8NF [21] for the NFO cross sections. These .K'-factors are 1.9 ± 0.2 for the two-photon process 
and 3.3 ± 0.5 for the three-photon process. The four-photon process is not included in NFO generators, 
and since the four-photon background is ~10~^ of the total background expectation in the inclusive signal 
region, the three-photon A'-factor is applied to the four-photon background sample. The uncertainties for 
these .K'-factors are determined by multiplying the renormalisation and factorisation scales independently 
by 2 and 0.5 and taking the largest deviations from the nominal value of the .fir-factor. 

The reducible backgrounds where electrons are misidentified as photons originate from multiple sources. 
Processes where a Z boson decays to an e'*'e~ pair, accompanied by a photon not from the matrix element, 
are modelled with Powheg-Box 1.0 [22], using Pythia 8 for fragmentation and hadronisation, and Z+y 
production is modelled with Sherpa 1.4.1 [23]. Backgrounds from processes involving the leptonic decay 
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of the W boson in association with photons and/or hadronic jets are simulated with Alpgen [24] and 
Herwig [25,26]. Possible mis-measurement of the rate of electrons misidentified as photons in simulation 
is addressed by comparing to electrons misidentified as photons from Z —> e'''e~ events in data. 

To obtain estimates of the rates at which true photons populate the regions of kinematic phase space 
assumed to be dominated by jets (used in the calculation of the systematic uncertainty for the data-driven 
estimate of jet backgrounds), a sample containing events with one hard-process quark or gluon and one 
prompt photon is simulated using Pythia and the CTEQ6L1 PDF set. 

The PDF sets for the simulated event samples of background processes used for the final background 
esfimafe in fhe inclusive search, for fhe MadGraph, Pythia and Alpgen -i- Herwig samples, are faken 
from CTEQ6F1, while for fhe Powheg-Box and Sherpa samples fhe PDF sefs are faken from CTIO. 


4.2 Simulated signal processes 

The Zyyy effecfive verfex has been implemenfed wifh FeynRules [27,28] and fhen used in a cusfomised 
MadGraph 5 model which is employed fo simulafe evenfs, using fhe CTEQ6F1 PDF sef and Pythia 8 
for fragmenfafion and hadronisafion. Each of fhe fwo non-frivial, independenf, lowesf-order effecfive 
Fagrangians for fhis process [29] confains a dimensionful coupling consfanf, and fhe values of fhese 
consfanfs have been calculated [30] using fhe SM expected Z —> 3y branching rafio of 5.41x10“^°. These 
SM values are used in fhe simulafion. The BSM process of a Higgs boson produced via gluon fusion and 
decaying fo four phofons via a pair of infermediafe a particles is simulafed wifh Powheg-Box and Pythia 
8 (using fhe CTIO NFO PDF sef). The BSM process of a new vector gauge boson decaying fo fhree 
phofons via Z' ^ a -i- y ^ 3y is simulafed wifh Pythia 8 (using fhe MSTW2008FO [21] leading-order 
PDF sef). 


4.3 Minimum-bias interactions and the ATLAS detector simulation 

Minimum-bias profon-profon inferacfions in fhe same or nearby bunch crossings (pileup) are modelled 
wifh Pythia 8, using fhe MSTW2008FO PDF sef. These pileup evenfs are overlaid onto fhe hard- 
scaffering process for all simulafed signal and background samples to reproduce fhe disfribufion of fhe 
average number of inferacfions per bunch crossing observed over fhe course of dafa-faking in 2012. 

All signal and background samples are processed wifh fhe full ATFAS defector simulafion [31] based on 
Geant 4 [32] and reconsfrucfed using fhe same soffware as fhaf used for collision dafa. 


5 Background composition estimate 

The backgrounds in fhe search for excesses in fhe inclusive signal regions are esfimafed from a com- 
binafion of simulated samples (defailed in fhe previous section) and mefhods employing collision dafa. 
The dominanf backgrounds in fhe inclusive signal region are fhe irreducible SM fwo-, fhree- and four- 
phofon processes, while for fhe Z —> 3y search channel, backgrounds involving elecfrons misidentified as 
phofons are dominant. 
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5.1 Backgrounds estimated from simulation 


The irreducible SM two-, three- and four-photon backgrounds, as well as backgrounds from processes 
involving electrons in the final state originating from Z decays and those involving the leptonic decay of 
the W boson in association with photons and/or hadronic jets, are estimated via simulation. The third 
photon for the SM two-photon background process typically arises from pileup interactions, but can 
occasionally be a quark- or gluon-initiated jet radiated from the incoming partons which is misidentified 
as a photon. Possible double-counting with the 2y + 1-jet final state (estimated via a data-driven method 
described in the following section) is avoided by omitting from consideration events in the SM two-photon 
simulated sample where one of the three photon candidates is a jet, using generator-level information. 
Possible mis-measurement of the rate of electrons misidentified as photons in simulation is addressed by 
comparing Z ^ e'*'e~ processes in simulation and in data. The per-electron scale factor is the ratio of 
the misidentification rate determined in data to that determined in simulated samples. This scale factor is 
independent of px and rj for the ranges considered here, and is found to be 1.03 ± 0.04. 


5.2 Data-driven estimates of ly + 1-jet, I7 + 2-jet, and 3-jet backgrounds 

A crucial aspect of the analysis is the data-driven estimate of the backgrounds where hadronic jets are 
misidentified as photons (hereafter called “jet fakes"), i.e., SM processes that can produce ly + 1-jet, I 7 
-I- 2-jet, and 3-jet events. Collision data are used to derive efficiencies for photons passing the isolation 
criterion (e^) and rates at which jets are misidentified as isolated photons (/jet)- These values of 6 y and ^et 
are then used in a likelihood matrix method (described below) to estimate the jet backgrounds. 

A sample of photon candidates consisting mainly of jet fakes is defined in fhe following way. The standard 
tight and loose photon identification categories are augmented with a medium definition [33], intermediate 
between tight and loose. The medium photon is defined by relaxing some EM shower shape requirements 
that provide high levels of rejection of jet fakes. When the medium definition is combined with a further 
requirement that the photon candidates fail tight (the combination hereafter called non-tight), the result is 
a sample of photon candidates that is primarily composed of jet fakes. This method presupposes that the 
distribution of the non-tight sample is composed primarily of jet fakes, and that the subset of tight 
photons with higher values of (the “tail", here for El^'^ > 1 GeV) is dominated by jet fakes. Under 
these assumptions, the tail of the non-tight distribution is scaled to match that of the tight distribution, 
thus providing a determination of the contribution of jet fakes to the signal region, i.e., the collection 
of photons that pass both the tight and the isolation criteria. The scaled non-tight distribution is then 
subtracted from the tight distribution. Photon isolation efficiency, €y, is then calculated as the ratio of 
the number of isolated photons (those that satisfy E}^° < 4 GeV) to the total number of photons in the 
tight distribution after this subtraction has been performed. The rate at which jet fakes are identified as 
photons, ^et^ is the ratio of the number of isolated photons to the total number of photons in the non-tight 
distribution. 

The assumptions described above are validated using simulated samples of events containing photons 
and jets, described in Section 4. Any collection of photon candidates consists of some combination of 
actual photons, which can be defined as “true", and other objects that are misidentified as photons. The 
non-zero true photon contamination in the set of non-tight photons, and the set of tight photons that fail 
the isolation criterion, is taken from the simulated samples and is used to derive a systematic uncertainty 
(described in Section 6 ) on the jet background estimate procedure. 
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Kinematic region Fraction satisfying isolation criterion 

Photons (e-y) Jets misidentified as photons (^et) 

1: 15 GeV < pj < 40 GeV, \rj\ < 1.37 0.939 ± 0.007 + 0.009 0.424 ± 0.001 ± 0.013 

2: pT >40GeV, It/I < 1.37 0.906 ± 0.006 + 0.013 0.256 ± 0.002 ± 0.010 

3: 1.52 < |?7| < 2.37 0.933 ± 0.007 ± 0.009 0.431 ± 0.002 ± 0.013 


Table 1; Photon isolation efficiencies (ey) and rates of jets misidentified as photons (^et) from collision data for the 
three kinematic regions, used for the jet background estimate. The isolation criterion is Ey° < 4 GeV. The three 
regions were chosen to maintain a large number of events in each bin. The first uncertainty is statistical while the 
second is systematic. 


The procedure is performed separately for three kinematic regions as follows. Photons are ordered by 
Pt. highest to lowest. Three regions in the px-p plane are defined as 1) 15 GeV < px < 40 GeV and 
Ipl < 1.37, 2) Px > 40 GeV and |p| < 1.37, and 3) 1.52 < |p| < 2.37. The separation into lower and higher 
Px bins around 40 GeV is chosen because this is the value at which ey and /jet are changing rapidly, and 
the three regions were chosen to maintain a large number of events in each bin. The values of ey and /jet 
are then calculated for each of the three regions, and the results are shown in Table 1 . 

The data-derived ey and /jet values are applied to events with three photon candidates to estimate the 
SM 27 + 1-jet, I 7 + 2-jet, and 3-jet backgrounds. This is done using a likelihood-based version of a 
standard matrix method (here called the “likelihood matrix method"). In standard matrix methods [33], 
a matrix of efficiencies relates an observed event that falls into a particular event category (based on 
some discriminating variable or variables) to the true, unknown final states to which the event has the 
possibility of corresponding, and the matrix is inverted to determine probabilities that a given observed 
event corresponds to one of these true final states. When summed over a large number of events, these 
per-event estimators average to the overall estimate of the number of events in each true final sfate. 

In the likelihood matrix method, by contrast (and with respect to the present three-photon search), the ex¬ 
pected yield for each three-object final state consisting of jets plus photons or all jets is the result of fitting 
a likelihood function to data. For the event sample where all three photons, ordered from highest to lowest 
Px, have satisfied fhe tight requirements, events are placed into 160 orthogonal categories designated by 
six criteria. These are defined first by the three regions in the px-p plane to which each photon candidate 
belongs. These are the same three regions that are used to categorise photons and to calculate ey and /jet, 
described and labeled previously as regions 1, 2, and 3. The remaining three criteria by which each event 
is categorized are three boolean variables, one for each photon candidate, indicating whether it passed or 
failed the isolation criterion. Since each of the three photon candidates either passes (P) or fails (F) isol¬ 
ation, there are 2^ = 8 possible isolation combinations for three photons: PPP, PPF, PFP, FPP, PFF, FPF, 
FFP, and FFF The three photons are ordered by px, from highest to lowest, and, since one of the three 
kinematic regions defined above depends only upon p, there are twenty possible px-p bin combinations: 
333, 332, 323, 322, 331,313,311,321,222,223,232, 233,221,211,231,213, 111, 113, 131, 133. This 
results in 8x20 = 160 categories, denoted PPP_333 for those events where the three photon candidates 
all passed isolation and had px-p values placing them in the “3" kinematic region, PPF_321 for those 
events where the leading and subleading photons passed isolation and the sub-subleading photon failed 
isolation, and the px-p value combinations placed them successively in the “3", “2", and “1" regions, 
etc. 
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Each of the 160 categories corresponds to a Poisson function where the observed number of events is 
the number of events seen in data for that category and the expected number of events is a sum of terms 
corresponding to each of the possible true (unknown) final states consisting of photons and jets or only 
jets for a particular pi-J] combination. Each term in a given sum is multiplied by the appropriate values 
of 6j, and ^et- A likelihood is then constructed consisting of a product of the 160 Poisson functions. The 
expectations for each true final slate are Ihe maximum likelihood esfimalors lhal resull from filling Ihis 
likelihood funclion fo fhe dala. Thai is, fhe frue unknown expecfafions are allowed lo floaf in Ihe fil and 
are consfrained lo be posilive and, hence, physical. The eslimaled number of evenls of a given final slate 
in a particular signal or conlrol region — defined by whelher Ihe photons passed or failed isolation — is 
determined by summing Ihe resulting expeclalions from Ihe fil limes Ihe appropriate €y and ^et values. 


6 Systematic uncertainties 

6.1 Data-driven uncertainties 

Eor the data-driven jet background estimate, systematic uncertainties arise in the calculation of the rate 
of photons passing the isolation criterion, ey, and the rate of jets misidentified as isolated photons (“jet 
fakes"), ^et- This calculation relies upon the assumptions that both the tail {E^^° > 1 GeV) of the E™ 
distribution of tight photons and the entirety of the E™ distribution of non-tight photons are primarily 
composed of jet fakes. Tests on simulations of photons and jets indicate that the true photon contamination 
in these jet-dominated regions is between 5% and 15%, depending on the region. These values are used to 
calculate different values of €y and /jet (where the number of photon candidates in a given region is altered 
by the corresponding percentage) which are then used in the jet background estimate. The deviations from 
the nominal signal region yield (assuming no true photon contamination) are calculated separately for the 
three final states of 2y -i- 1 jet, ly -i- 2 jets, and 3 jets, and these values (4%, 10% and 21%, respectively) 
are taken as systematic uncertainties on the estimates of these backgrounds in the signal region. 

An additional uncertainty associated with the data-driven methods employed arises from the choice of 
kinematic variable used to categorise photons and then to calculate and apply ey and ^et- The baseline 
analysis uses three bins in the pj-p plane, described in Section 5, and separate analyses using either px- 
dependence only or p-dependence only are conducted as well. The largest deviation of the two different 
methods from the nominal method is 13%, which is taken as a systematic uncertainty. 


6.2 Simulation uncertainties 

The uncertainty on the integrated luminosity for the data sample is 2.8%, derived using the same meth¬ 
odology as that detailed in Ref [34]. 

The photon identification efficiency has been direcfly measured in dafa using photons from Z ^ e'^e~ /p'^p 
radiafive decays [12]. The sysfemafic uncerfainfies on fhe signal region yield due to fhe uncertainly on 
Ihis efficiency measuremenl are found lo range from <1% fo 6% for simulated backgrounds, and from 
3% lo 7% for simulaled signal processes, depending on the sample. 

As mentioned in Section 3, the analysis supplements the isolation prescription — Ej° < 4 GeV, with a 
cone size of AR < 0.4 — with an isolation energy correction that is applied to photons with overlapping 
isolation cones. This procedure improves sensitivity to lower-mass two-photon resonances where the 


photon pairs are close together in AR. To account for possible over- or under-correction due to a photon 
being near the edge of the isolation annulus, an additional systematic uncertainty is assessed. The pj 
values of all isolated photons in simulated samples are calibrated to yield agreement with the values 
observed in data [35]. Since the calibration factors for isolated photons deviate from one by typically 
less than 5%, a value of 5% is a conservative estimate of the uncertainty on photon pj. To assess the 
systematic uncertainty on the isolation energy correction, the measured value of the p^ of the other tight 
photon in the isolation cone is varied by 5%, the correction procedure is applied, and the effects are 
propagated to the final event selection in the signal region. For example, using simulated samples of a 
Higgs boson decaying to four photons via H ^ aa ^ 4y, the systematic uncertainty due to this effect is 
smaller for higher ratios of nialniH (as large as 6% when the pT is varied by -5% and <1% when the pT is 
varied by +5%, for mn = 900 GeV and nia = 440 GeV), and the uncertainty is larger for smaller ratios of 
malniH (as large as 69% when the pj is varied by -5% and 12% when the px is varied by -i-5%, for niH = 
900 GeV and nia = 50 GeV), as the photons tend to overlap within the isolation cone more frequently. 

The uncertainties on the event yields due to systematic uncertainties in the photon energy scale and 
resolution [35] are found to range from <1% to 4% for the simulated signal and background samples. 
The uncertainties on the event yields due to systematic uncertainties in the scale factors used to yield 
agreement between photon identification efficiencies calculated in data and simulated samples [35] are 
found to range from <1% to 8% for simulated backgrounds, and from 1% to 4% for simulated signal 
processes. The systematic uncertainty on the scaling factor for electrons misidentified as photons in 
simulated samples is taken to be the statistical uncertainty arising from the calculation, i.e., 4%, since this 
is as large as or larger than the systematic uncertainties due to the photon energy scale and resolution, 
above. The efficiency and uncertainties of the three-photon trigger chain have been determined to be 
98.5% ± 0.1% (stat.) ± 0.2% (syst.). The trigger efficiency is cafcufated using singfe photons (with px 
vafues corresponding to the vafues used for the analysis event selection) from Z boson radiative decays 
and then, under the assumption that the per-event performance of the photon trigger for one photon is 
uncorrelated to that for another photon in the same event, multiplying these values to obtain the overall 
trigger efficiency. 

Uncertainties on calcufated cross sections for simulated background processes due to QCD renormal¬ 
isation and factorisation scales and due to the choice of PDF set and value of as used in simulation are 
addressed via the recommendations of PDF4LHC [36]. The resulting combined uncertainties are found to 
range from 7% to 16%, depending on the simulated background process. The total theoretical uncertainty 
on the SM 3y background process due to the uncertainty on the LO to NLO correction, combined with 
the uncertainties due to choice of PDF set and renormalisation and factorisation scales, is found to be 
-30%. 

Uncertainties exist for the measured or calculated production cross sections for SM particles for which 
BSM decays are considered as signal scenarios and are accounted for. For the BSM Higgs boson scen¬ 
ario of h ^ aa ^ 4y the gluon fusion production cross section for the SM Higgs boson with m/, = 
125 GeV [37,38], cr/,,sM - 19.27 pb is used, with an uncertainty of ±10.4% due to choice of PDF set and 
renormalisation and factorisation scales. For the rare decay Z ^ 3y, the measured pp ^ Z cross section 
of (2.79 ± 0.02 ± 0.11) X 10“^ pb [39] is used. An additional uncertainty of ±12.3% — determined by 
varying the QCD renormalisation and factorisation scales, PDF set, and value of as — is also assessed, 
to account for variations in the simulation of the kinematics of the final-state photons and, hence, the 
acceptance in the signal region. 

These systematic uncertainties are summarized in Table 2. 
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Dafa-driven 

Background 

Signal 

Phofon confaminafion of confrol regions 

4-21 

— 

Kinematic paramefrizafion 

13 

— 

Simulation 

Phofon ID 

1-6 

3-7 

Phofon isolation correcfion 

<1^ 

<1-69 

Phofon energy scale and resolution 

<1^ 

<1^ 

Phofon scale faclors 

<1-8 

1-4 

Elecfron scale faclors 

4 

4 

Trigger 

0.2 

0.2 

Luminosily 

2.8 

2.8 

Cross section 

7-16 

4-12 


Table 2: Systematic uncertainties (%) on the expected event yields in the signal region. The values given for data- 
driven backgrounds correspond to the three jet backgrounds described in the text. For simulated samples, when a 
range is given it corresponds to the smallest and largest uncertainties for all simulated backgrounds or signals. 


7 Results and interpretations 

This section presents results based on the number of events in a broad inclusive region and a restricted 
region focusing on the rare decay Z —> 3y, as well as results from the search for resonances in the 
di-photon and tri-photon invariant mass spectra. 


7.1 Inclusive and Z 3y regions 

The number of SM background events expected in the signal region is 1370 ±140 (combined statistical 
and systematic uncertainties) and the observed number of events is 1290. The observation is in agreement 
with the SM expectation. Additionally, while the event selection is optimised for a search for physics 
beyond the SM as opposed to a measurement of the 3y inclusive cross section, the results are nevertheless 
compatible with the irreducible all-photon process expectations from the SM. The expected and observed 
yields in the signal region are presented in Table 3, and the expected and observed yields in signal and 
control regions where all three photons have passed the tight identification criterion are shown in Figure 2. 
In the figure, fhe red hafched band, in fhe signal region bin, is fhe combination of sfafislical and sysfemafic 
uncerfainfies on all background sources, while fhe black hafched band, in fhe confrol regions, is fhe 
combination of sfafisfical uncerfainfies of fhe dafa-driven jef background esfimafe and fhe expecfed yields 
from simulafed samples of SM background processes. For fhe inclusive signal region, fhis corresponds 
fo a model-independenf observed (expecfed ±lcr) upper limif, af fhe 95% confidence level (C.L.), on 
fhe number of signal evenfs of 240 (273^^^), and fo fhe model-dependenf upper limifs on fhe inclusive 
fiducial cross secfion in fhe aforementioned accepfance for fhe signal scenarios of fhe BSM Higgs boson 
and Higgs boson-like decays and fhe Z' decays shown in Table 4, where hypofhesis fesfing and limit 
setting are calculated using the profile likelihood ratio as fhe fesf sfafisfic for fhe CLs fechnique [40] in 
fhe asympfofic approximation [41,42]. The fiducial efficiencies for each signal scenario are defermined 
wifh respecf fo a generafor-level kinematic region wifh fhe same requiremenfs applied fo fhree-phofon 
evenfs as fhose used for fhe analysis signal region. This fiducial region is defined as fhe sef of evenfs fhaf 
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Figure 2: Observed and expected yields in signal and control regions for the full mass range (left) and the restricted 
range of 80 GeV < OT 3 y < 100 GeV (right), for events where all three photon candidates satisfy the tight photon 
identihcation criteria. The bins along the horizontal axis correspond to orthogonal subsets of events where each 
subset is categorised by whether the three photons — ordered from largest to smallest values of px— passed (“P") or 
failed (“F") the isolation criterion. The leftmost bin is the signal region, composed of events satisfying PPP, and the 
other bins are the different control regions, where at least one of the photon candidates failed the isolation criterion. 
The red hatched band, in the signal region bin, is the combination of statistical and systematic uncertainties, while 
the black hatched bands represent statistical uncertainties. As a result of the data-driven jet background estimate, 
the statistical uncertainty in each bin is partially correlated with the uncertainty on the data in that bin. 


contain three photons where 1) each photon satisfies a pseudorapidity requirement of I 77 I < 2.37, excluding 
the transition region between the barrel and end-cap of 1.37 < \r]\ < 1.52, 2) the three photons satisfy pj > 
22, 22, and 17 GeV, and 3) each photon satisfies < 4 GeV, where is a generator-level 

definition of the photon isolation criterion equivalent to that used for event selection on reconstructed 
events. The fiducial efficiencies are similar for the considered signal scenarios for mass points where 
the distributions of photon pT^ for all photons, tend to peak higher than pj > 50 GeV. This is because 
the overall photon identification efficiency decreases for photons with pj < 50 GeV [12]. Since the pj 
distribution for at least one of the photons for signal scenarios with lower-mass resonances tends to peak 
at lower values, the fiducial efficiencies are lower. 

Using the same data-driven and simulation-based methodology restricted to the region 80 GeV < ^ 3 ^ < 
100 GeV provides a test for the rare decay of the Z boson to three photons. The SM branching ratio for 
the process is predicted to be ~5xl0~^° [9], but it has yet to be observed. Table 3 (right) and Figure 2 
(right) summarise the observed counts as well as background expectations in this restricted region. The 
data are consistent with the SM expectation: 244 events are observed and 233 ± 28 events are expected, 
while the signal expectation from simulation, for BR(Z ^ 37 ) = 10“^ (corresponding to the previous 
limit from LEP [10]), is 418 ± 9 events. Using the same hypothesis-testing and limit-setting procedure 
described above, but taking the signal expectation from the simulated sample described in Section 4, the 
observed (expected) limit, at the 95% C.L., on the branching ratio of the Z boson decay to three photons 
is found to be BR(Z —> 3y) < 2.2 (2.0) xl0~^, a result five times stronger than the previous LEP limit of 
< 10 -^ 
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Process 

Inclusive signal region 

80 GeV < m^y 

< 100 GeV 

2y (sim.) 

330 ± 

50 

24 ± 

8 

3y (sim.) 

340 ± 

110 

30 ± 

10 

4y (sim.) 

1.3 ± 

0.4 

0.07 ± 

0.02 

2y,lj (D-D) 

350 ± 

60 

65 ± 

19 

lr,2j (D-D) 

110 ± 

40 

13 ± 

10 

3j (D-D) 

43 ± 

11 

6.1 ± 

2.0 

Z —> e'''e~ (sim.) 

+1 

CO 

22 

43 ± 

13 

Z+y (sim.) 

89 ± 

11 

+1 

CO 

6 

VF-i- 7 -i-( 0 ,l, 2 )j (sim.) 

11.4± 

1.5 

2.7 ± 

0.7 

VF-i-2'y-i-(0,l,2)j (sim.) 

6.1 ± 

0.5 

0.68 ± 

0.08 

Tolal SM exp. 

1370 ± 

140 

233 ± 

28 

Observed 

1290 


244 



Table 3; Expected and observed event yields in the inclusive signal region and for the signal region with a further 
requirement of 80 GeV < m^y <100 GeV. Background expectations estimated via simulations are marked sim., 
whereas data-driven calculations are denoted as D-D. The uncertainties for each row are the combination of statist¬ 
ical and systematic uncertainties for a given background process, and the overall uncertainties in the second to last 
row are the combined uncertainties for the total background expectations for each signal region. 


7.2 The ly and 3/ resonance searches 

In addition to the tests based on the number of events in the inclusive signal regions, searches are per¬ 
formed for resonances in the two-photon and three-photon invariant mass {m 2 y and m^y) distributions 
for events in the inclusive signal region. For these resonance searches, the background contribution is 
estimated from a fit to the m 2 y or m^y sideband regions, and thus does not rely upon simulated samples 
for the background estimate. The sideband is modelled as a fourth-order polynomial, and the size of 
the sideband is mass-dependent, symmetric around the hypothesised resonance mass, following a local- 
spectrum approach. The range of the observed mass spectrum that is used for the sideband fit is a local, 
truncated subset of the full spectrum. For the m 2 y {mT,y) resonance search, the sideband is 20 (25) GeV in 
each direction for m 2 y (mjy) < 90 (230) GeV, where the event counts change rapidly as a function of m 2 y 
(m^y), and rises to a sideband size of 80 (100) GeV in each direction for m 2 y {mi,y) > 195 (425) GeV, in¬ 
creasing roughly linearly with mass as the spectrum becomes smoother. The m 2 y (m^y) resonance search 
begins at a mass hypothesis of 10 (100) GeV, and proceeds in steps of 0.5 GeV. The signal component 
of the resonance search is a Gaussian function with a fixed widfh fhaf varies wifh particle mass, and fhe 
widfhs are defermined from simulafed signal samples. Since fhe simulafed signal samples are generafed 
wifh a narrow-widfh approximation for bofh fhe pseudoscalar a and fhe Z' in all cases, fhe 2y and 3y 
mass resolufions for fhis search are equivalenf fo Gaussian funcfions fhaf accounf for defecfor resolu- 
fion, and are defermined via fifs fo fhe simulafed signal samples. Hypofhesis fesfing and limif selling are 
performed using fhe profile likelihood ralio as fhe lesl slalisfic for fhe CL^ fechnique in fhe asymplofic 
approximalion. 

The resonance search is performed separalely for fhe Ihree Iwo-pholon mass speclra defined by fhe Ihree 
possible phofon pairings for Ihree photons in fhe inclusive signal region. As menfioned previously, fhe 
pholons are ordered by pj, from highesl fo lowesl, and so fhe Ihree Iwo-pholon mass speclra are denofed 
w^l 2 ^ w^i 3 and m 23 , where fhe 1, 2, and 3 refer fo fhe px-ordered pholons. The observed m 2 y and m^y 
speclra in fhe inclusive signal region are shown in Figure 3. Also shown in Figure 3, for visualisation 
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efficiency 

CTfid X 54 [fb] CToverall [fb] 
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nihiH [GeV] ma [GeV] 
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62 
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{ 21 * 1 ) 

1 118(155!?3) 

300 
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2 i:f) 

1 29(35!^) 
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100 

0.688 + 0.003 

17 1 

20-6) 

1 27(34+^) 
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100 

0.680 + 0.003 

17 1 

[2011) 

1 27 (33!^ 


Z' —> a + y —> 3y 
mz' [GeV] [GeV] 


100 

40 

0.438 + 0.009 

27 ( 

31!?) 

1 316 ( 

;387!?f) 

200 

100 

0.611 ±0.005 

19 1 

223 ) 

1 53 ( 

:62!?^) 

1 

400 

100 

0.649 ± 0.004 

18 ( 

717) 

1 51 ( 


1 

600 

100 

0.667 ± 0.004 

18 ( 

20 !^) 

1 39 ( 

;48+j2) 

1 

1000 

100 

0.636 ± 0.004 

19 1 

, 21 !^) 

1 38 ( 
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Table 4; Top row: Observed and expected model-independent upper limits on event yields for new physics processes 
for the inclusive signal region. Also shown are the efficiencies for the fiducial kinematic region defined in the text for 
some example mass points for the signal scenarios explicitly considered here, and the corresponding observed and 
expected (+lcr) upper limits on the fiducial cross section within the acceptance. Total statistical uncertainties are 
quoted for the fiducial efficiencies, and the uncertainties for the upper limits correspond to the uncertainties arising 
from the +lcr upper limits calculated via hypothesis testing using the combination of statistical and systematic 
uncertainties. 


purposes only, is the background expectation per bin, determined from the sideband fit to data as a part 
of the resonance search. The resonance search is performed with a step size of 0.5 GeV and so the final 
results shown in Figures 4 and 5 demonstrate sensitivity to resonances with widths appropriate to the BSM 
models considered here. The widths of the bins in Figure 3 do not correspond to the mass resolution for 
the signal scenarios in question. The background estimates and significances shown in Figure 3 provide 
a complementary comparison of the local agreement between data and expectation. The lower panels 
show the significance, in units of standard deviations of a Gaussian function, of the observation in each 
bin, taking into account the fractional uncertainty on the background as a result of the sideband fit. The 
significances shown in the lower panels in Figure 3 are derived from the p-value for the background- 
only hypothesis for each bin, calculated using a frequentist binomial parameter test [43-45]. For regions 
beyond the sensitivity of the search, no background estimate is shown. 

For the Hjh aa ^ Ay BSM signal scenario, the photon pairing (among the three pj-ordered photons) 
that most often corresponds to the photons arising from the decay of the same a particle is (2,3). As a 
result, the resonance search in the ^23 spectrum provides the best sensitivity to this model. The widths of 
the Gaussian signal component — corresponding to the detector resolution — are taken from simulated 
samples of a Higgs boson decaying to four photons via a pair of intermediate pseudoscalar a particles, 
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and vary from 0.6 GeV < croauss < 3.2 GeV for 10 GeV <ma < 440 GeV, and are largely independent of 
mh/H- 

No excess above background is detected, and upper limits, for a SM-like Higgs boson of m/, = 125 GeV 
(and assuming kinematics associated only with gluon fusion SM Higgs boson production), are calculated. 
Additionally, limits are set for Higgs boson-like scalars with masses larger than 125 GeV. The results of 
these resonance searches are shown in Figure 4 for the SM-like Higgs boson of nth = 125 GeV (in the top 
row) and, as an example of a higher scalar mass, for mh = 600 GeV (in the bottom row). The resonance 
search limits for higher values of mh are limited by the small number of events in the mass spectra 
at higher values of m 2 y As shown in Figure 4, the limits vary as a function of two-photon invariant 
mass, but an overall upper bound on the limits is determined to be cr x BR(/j ^ aa) x BR(a ^ yy)^ 
< 1 X 10“^a■sM^ for 10 GeV < nia < 61 GeV for the SM-like Higgs boson of nth = 125 GeV and, for 
the higher scalar mass case, cth x BR(// ^ aa) x BR(a —> yy)^ < 0.02 pb for lower nia values in the 
range 10 GeV < nia <90 GeV and <0.001 pb for higher nia (up to 245 GeV for the resonance search in 
the m 23 spectrum for inn - 600 GeV shown in Figure 4). Additionally, using the expected signal yields 
from simulated samples, inclusive limits are calculated for 300 GeV < niH < 900 GeV and for a range of 
nia < mnll, including values beyond the range of the mass spectra used for the resonance search. These 
inclusive limits are shown in Table 5. 

Moreover, upper limits on the Z' production cross section times the product of branching ratios, crz' x 
BR(Z' —> a -I- y) X BR(a —> yy), are found to be in the range of 0.04 pb to 0.3 pb, depending upon mz' 
and nia- Upper limits, at the 95% C.L., on crz' x BR(Z' —> a -i- y) x BR(a —> yy) are shown in Table 6, as 
a function of nia, using the expected signal yields from simulated samples. Additionally, using a narrow- 
width approximation to the Z' resonance width, local excesses corresponding to Gaussian resonances due 
to detector resolution are searched for in the m^y spectrum. The Gaussian widths are determined via fits 
to the Z' simulated signal samples, described in Section 4. For the range of my for which the resonance 
search is possible, the Z' width exhibits a small dependence on For each Z' mass point, three different 
samples are simulated, with different values of ma- The average of the three measured Z' widths for each 
of the ma points simulated is taken as the width for a given my , and these values are used for the three- 
photon resonance search, interpolating for niy points between those for which samples are simulated. 
These values range from 1.5 GeV < crcauss < 2.4 GeV for 100 GeV < niz’ < 500 GeV. The results, along 
with the local p-values for the background-only hypothesis, are shown in Figure 5. The smallest local 
p-value is found to be 0.0003 (3.4cr local significance), at mj,y = 212 GeV which, after adjusting for a 
trials factor [46], corresponds to a global p-value of 0.087 (1.4cr global significance). 
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Figure 3: Observed spectra of m\ 2 , '«i 3 , and OT 23 , where the 1, 2, and 3 refer to three pj-ordered photons, as well 
as m^y. For illustration purposes only, also shown is the expected background per bin, determined via unbinned 
sideband fits to the data as a part of the resonance search, for a hypothesised resonance mass defined by the centre 
of the bin, as well as the signal expectation for a few mass points for the BSM scenarios considered here. The lower 
panels show the significance, in units of standard deviations of a Gaussian function, of the observation in each bin, 
taking into account the fractional uncertainty on the background as a result of the sideband fit. This significance is 
derived from the p-value for the background-only hypothesis for each bin, calculated using a frequentist binomial 
parameter test [43-45]. The signal distributions used for the m 2 y resonance searches have two components, a 
narrow Gaussian core for correctly paired two-photon combinations and a wide distribution for incorrectly paired 
combinations that is well described by the polynomial used to simultaneously model the background shape for the 
resonance search described in Section 7.2. 
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Observed (expeeted) 95% C.L. upper limits 

niH [GeV] 

nta [GeV] 


(Th X BR(// ^ aa) x BR(a ^ yy)^ [fb] 


20 

50 

100 140 

300 

48(60:|5) 

33(40!|) 

29 ( 352 ^) 28 (342^) 


50 

100 

200 290 

600 

31 (38+^0 j 

27(34+7) 

25(312^6) 25(312^) 


50 

100 

200 440 

900 

36(44tV) 

27 (3327) 

26 (332^6) 26 (3227) 


Table 5: Expected and observed 95% C.L. upper limits on cth x BR(H —> aa) x BR(a —» The uncertainties 

for the expected limits are the +lcr uncertainties resulting from the hypothesis tests for each mass point, taking into 
account statistical and systematic uncertainties. 


Observed (expeeted) 95% C.L. upper limits 

mz' [GeV] 

nia [GeV] 



(Tz' X BR(Z^ 

^ a + y) X BR(a 

^ yy) [fb] 


40 

60 

80 

100 

320(3902f§) 

150 (1702|[]) 310 (370+_i80oO) 


50 

100 
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200 

78 (902i) 

53(62«») 

51 ( 58 : 13 ) 
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51 (6321^) 

44(55+_\2) 

38 (47212) 
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39(48+_i92) 

4l(52+_io) 

4l(52+_V) 


100 

350 

600 

800 

38 (46+_V) 

35 (4321) 

35(4321) 


100 

450 

800 

1000 

38 (46+_V) 

54(6421^) 

37 (43+_i82) 


Table 6: Expected and observed 95% C.L. upper limits on cr^' x BR(Z' —> a + y) x BR(fl —> yy). The uncertainties 
for the expected limits are the +lcr uncertainties resulting from the hypothesis tests for each mass point, taking into 
account statistical and systematic uncertainties. 
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Figure 4; Left; Local p-values for the background-only hypothesis as a result of a resonance search with respect to 
the BSM process hjH ^ aa ^ 4y, for m/, = 125 GeV (top row) and mn - 600 GeV (bottom row), as a function 
of nia, determined via a search for local excesses in the OT 23 spectrum. Right; Upper limits, at the 95% C.L., on 
(cr/crsM)xBR(/! ^ afl)xBR(fl —> yy)^ (top row) and ct/^xBR(// —> fla)xBR(fl —> yy)^ (bottom row). Also shown 
are the +1 and 2 cr uncertainty bands resulting from the resonance search hypothesis tests, taking into account the 
statistical and systematic uncertainties from simulated signal samples which are used to determine signal efficiency 
and Gaussian resonance width due to detector resolution for each mass hypothesis. 
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Figure 5; Left; Local p-values for the background-only hypothesis as a result of a resonance search with respect 
to the production of a new vector gauge boson Z' as a function of niz', determined via a search for local excesses 
in the m-iy spectrum, using a narrow-width approximation to the Z' resonance width. The smallest local p-value 
is found to be 0.0003 (3.4cr) which corresponds to a global p-value of 0.087 (L4cr). Right: Upper limits, at 
the 95% C.L., on ctz' x BR(Z' —» a + y) x BR(fl ^ yy). Also shown are the +1 and 2cr uncertainty bands 
resulting from the resonance search hypothesis tests, taking into account the statistical and systematic uncertainties 
from simulated signal samples which are used to determine signal efficiency and Gaussian resonance width due to 
detector resolution for each mass hypothesis. 
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8 Conclusion 


A search for new phenomena in events with at least three photons has been performed using 20.3 fb“^ of 
LHC pp collision data at = 8 TeV collected with the ATLAS detector at CERN. The SM background 
expectation is in agreement with the data, and is determined to be 1370 ± 140 events while 1290 events 
are observed. The model-independent observed (expected) 95% C.L. upper limit on the number of signal 
events is found to be 240 (273^gg). Upper limits at the 95% C.L. are calculated on the fiducial cross 
section ctm for events from non-SM processes for several signal scenarios. The observed (expected) limit 
on the branching ratio of the Z boson decay to three photons is found to be BR(Z ^ 3y) < 2.2 (2.0) 
xl0“^, a result five times sfronger fhan fhe previous resulf from LER 

In addifion, a search for local excesses in fhe fwo-phofon and fhree-phofon invarianf mass disfribufions 
is conducted. Eor fhe fwo-phofon mass specfra, no significanf excesses are defecfed, and fhe 95% C.E. 
upper limif on (ct/ctsm) x BR(/i ^ aa) x BR(a ^ yy)^ (assuming kinemafics associated only wifh gluon 
fusion SM Higgs boson production) is calculated fo vary from ~3xl0“^ fo ~4xl0~"* for 10 GeV < nia < 
62 GeV for a SM-like Higgs boson wifh a mass of nih - 125 GeV. Eimifs are sef for Higgs boson-like 
scalars H wifh masses up fo m// = 900 GeV and are found fo be cth x BR(// ^ aa) x BR(a ^ yy)'^ < 
0.02 fo 0.001 pb, depending upon mn and m^. Eor fhe fhree-phofon mass specfrum, fhe resonance search 
is conducted in fhe confexf of a Z' decaying fo fhree phofons. The smallesf local p-value is found fo be 
0.0003 (3.4cr local significance), af mz' - 212 GeV which, afler adjusfing for a frials factor, corresponds 
fo a global p-value of 0.09 (L4cr global significance). Upper limifs af fhe 95% C.E. on fhe Z' production 
cross secfion times fhe producf of branching rafios, crz’ x BR(Z' a+ y)x BR(a ^ yy), are found fo be 
in fhe range of 0.04 pb fo 0.3 pb, depending upon mz’. 

These model-independenf resulfs are fhe firsf of fheir kind, as are fhe inferprefafions for a Higgs boson 
decaying fo four phofons via fwo intermediate pseudoscalar a particles (for a SM-like Higgs boson of 
m/, = 125 GeV and for Higgs-like scalars of higher masses) and for fhree-phofon resonances correspond¬ 
ing fo a new vecfor gauge boson. 
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Y Eang^^®, M. Eanti^^^’^'*’, A. Earbin^, A. Earilla^^^^, T. Earooque'^, S. Earrell^^, S.M. Earrington^™, 
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Y. Guo^^’’, S. Gupta^^*’, G. Gustavino^^^^d 32 b^ p Gutierrez'N.G. Gutierrez Ortiz’^, C. Gutsebow^^, 
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K. A. Johns^, K. Jon-And^^^^’^^^'’, G. Jones'^*^, R.W.L. Jones^^, T.J. Jones^^, J. Jongmanns^*^, 
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S. Kohlmann'^^, Z. Kohout'^^, T. Kohriki®, T. Koi'®, H. Kolanoski'^, I. Koletsou^, A.A. Komar®’*, 

Y. Komori'^^, T. Kondo®, N. Kondrashova®, K. Kbneke®, A.G. Kbnig'®, T. Kono®, 

R. Konoplieh""’", N. Konstantinidis^^, R. Kopeliansky'^^, S. Koperny^*^, E. Kbpke^^, A.K. Kopp®, 

K. Koreyl^^, K. Kordas'^^, A. Korn^*, A.A. Korol'®’^, I. Korolkov'^, E.V. Korolkova'^^, O. Kortner"", 

S. Kortner"", T. Kosek'^^, V.V. Kostyukhin^', V.M. Kotov®, A. Kot'wal®, 

A. Kourkoumeli-Charalampidi'^^, C. Kourkoumelis^, V. Kouskoura^^, A. Koutsman'^^^^, 

R. Kowalewski'®, T.Z. Kowalski^^®, W. Kozaneeki'^^, A.S. Kozhin'^", V.A. Kramarenko®, 

G. Kramberger^^, D. Krasnopevtsev^^, M.W. Krasny^", A. Krasznahorkay^", J.K. Kraus^', 

A. Kravehenko^^, S. Kreiss"", M. Kretz^^'^, J. Kretzsehmar^^, K. Kreutzfeldt^^, P. Krieger'^*, 

K. Krizka^', K. Kroeninger®, H. Kroha"", J. Kroll'^^, J. Kroseberg^', J. Krstie'^, U. Kruehonak^^, 

H. Kriiger^', N. Krumnaek®, A. Kruse'^^, M.C. Kruse®, M. Kruskal^^, T. Kubota*^, H. Kueuk’^, 

S. Kuday®, S. Kuehn^^, A. Kugel^^^^, E Kuger'^"', A. Kuhl'^^, T. Kuhl®, V. Kukhtin®, R. Kukla'^^, 

Y. Kulehitsky®, S. Kuleshov^^'’, M. Kuna'^^®’'^®, T. Kunigo^^, A. Kupeo'®, H. Kurashige®, 

YA. Kuroehkin^^, V. Kus'®, E.S. Kuwertz'®, M. Kuze'^^, J. Kvita"^, T. Kwan'®, 

D. Kyriazopoulos'^^, A. Ea Rosa'^^, J.E. Ea Rosa Navarro®", E. Ea Rotonda^’^’®'’, C. Eaeasta'®, 

E. Eaeava'^^‘‘’'^^'’, J. Eaeey^^, H. Eaeker'^, D. Eaeour^", V.R. Eaeuesta'®, E. Eadygin^^, R. Eafaye^, 

B. Eaforge*", T. Eagouri'^^, S. Eai^^, E. Eambourne^^, S. Eammers^', C.E. Eampen’, W. Eampl^, 
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K. Eantzseh^', A. Eanza'^'®, S. Eaplaee^", C. Eapoire^", J.E. Eaporte'^", T. Eari^'^*, 
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27 



C. Leggett^^, A. Lehan^^, G. Lehmann Miotto^*’, X. Lei^, W.A. Leight^^, A. Leisos^^^’“', A.G. Leister^^^, 
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E. Eobodzinska"^^, P. Eoeh^, W.S. Eoekman^^’, E.K. Eoebinger^"^, A.E. Eoevsehall-Jensen^^, 
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U. Mallik^^, D. Malon^, C. Malone^^^, S. Maltezos^°, V.M. Malyshev^*’^, S. Malyukov^°, J. Mamuzie^^, 
G. Maneini^^, B. Mandelli^*’, E. Mandelli^^®, 1. Mandic^^, R. Mandryseh^^, J. Maneira^^^^’^^^’’, 
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